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We consider Reissner-Nordstrom solutions in quasi-topological gravity, obtaining exact solutions 
to the field equations yielding charged quasi-topological black holes. We study their thermodynamic 
behaviour over a range of parameters that yield ghost-free and stable space times. We find that 
a sufficiently negative quasi-topological parameter can yield black holes with 2 horizons, even for 
zero charge. We discuss the thermodynamic stability for the class of solutions we obtain. We also 
describe the structure of exact charged solutions to k th order quasi-topological gravity. 

I. INTRODUCTION 

Gauge/gravity duality - the notion that gravitational dynamics in a given dimensionality is dual to some other 
(nongravitational) field theory of a lower dimensionality - has been a subject of fruitful investigation in recent years. 
Its chief instantiation has been via the AdS/CFT correspondence PQ, and investigations of this relationship from 
the perspective of Einstein gravity via the trace anomaly [2] have indicated that the duality between central charges 
and couplings on the non-gravitational side and the parameters on the gravitational side takes place only for those 
conformal field theories for which all the central charges are equal. Since this is because there are not enough free 
parameters in Einstein gravity, there has been motivation to broaden the universality class of dual field theories to 
more general theories of gravity containing more free parameters. Examples are Lovelock theory [3] and various forms 
of quasi-topological gravity [MB], whose additional central charges have recently been investigated holographically 
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We are concerned here with cubic quasi-topological gravity coupled to electromagnetism in D dimensions. Such 
higher-order curvature terms modify the gravitational action in a way that may exhibit useful dualities, in a manner 
similar to the AdS/CFT correspondence pQ. For example, requiring the dual gauge theory to respect causality (by 
considering graviton fluctuations in the gravity theory, which probe the bulk geometry in the presence of a black 
hole) yields constraints on the coupling constants of the gravity theory so that the dual theory does not support 
superluminal signals. In Lovelock gravity these causality constraints precisely match those arising from requiring 
positive energy fluxes [HI HD] ■ However, this matching does not appear in general, in particular for cases where the 
gravitational equations of motion are not second order |llj . Cubic quasi-topological gravity is of interest insofar 
as there are three constraints that arise from requiring positive energy fluxes, which determine the three coupling 
constants. No evidence for causality violation has been observed when the curvature-cubed couplings are consistent 
with these constraints [7] . 

Quasi-topological gravity is akin to Lovelock gravity insofar as the field equations are second order differential 
equations with respect to the metric, though it produces nontrivial modifications to the action in dimensions where 
Lovelock gravity is a topological invariant [4]. However, these additional dimensions come with the caveat that 
the field equations are second order only when the solution is spherically symmetric. The nomenclature "quasi- 
topological" is employed because in the dimensions for which its terms do not normally affect the action, breaking 
spherical symmetry causes it to contribute nontrivially to the action (and therefore, quasi-topological terms are not 
true topological invariants), yielding 3 rd order differential equations. Despite this, one of the more intriguing features 
of quasi-topological gravity is that the linearized equations describing gravitons propagating in an AdS vacuum match 
precisely the second-order equations of Einstein gravity. It is therefore possible to find stable vacua in the theory that 
are free of ghosts [5]. 

Here we seek charged Reissner-Nordstrom (RN) black hole solutions to the field equations of quasi-topological 
gravity, and analyze their behaviour. RN black holes are significant in Einsteinian gravity, featuring a characteristic 
dual event horizon because of their charge. As in the Lifshitz case [12], all the desirable properties of quasi-topological 
gravity are preserved upon coupling to a Maxwell field, and so it is reasonable to explore the physics of charged 
solutions. Indeed, our motivation for the study of these black holes is the utility of their charge parameter in 
the context of gauge-gravity duality. The goal in a search for dual theories is the development of a "dictionary" 
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of parameters, which allows one to evaluate quantities in one theory and map them to another. Work has been 
performed which hints towards dualities between higher curvature theories and exotic condensed matter behaviour, 
such as the proposed relationship between Gauss-Bonnet gravity and non- Fermi liquids [T3]- In addition, studies 
have been performed to examine relationships between uncharged quasi-topological gravity with probe-limit fields 
and holographic superconductors [14] . 

It would be particularly enticing to perform similar work in quasi-topological gravity for a number of reasons. First, 
the third-order quasi-topological terms become active (nontrivial) in five dimensions, allowing us to hypothesize a 
dual between a five dimensional quasi-topological theory and a four dimensional physical theory, with the bonus of a 
tuneable quasi-topological parameter. In addition, the seemingly unnatural constraint of spherical symmetry in quasi- 
topological gravity becomes unimportant, because the vast majority of metrics we understand well are spherically 
symmetric, so a search for spherically symmetric dual behaviour is likely to be equally profitable to a more general 
search. 

RN black holes feature a charge parameter in addition to their mass, which may be required to reproduce a 
constant low-temperature Sommerfeld ratio (as in Gauss-Bonnet gravity [13.). The Sommerfeld ratio's behaviour is a 
fundamental property of non-Fermi liquid systems, and recovering a constant Sommerfeld ratio at low temperature is 
a start to the development of a holographic duality. We will explore the quasi-topological behaviour of charged black 
holes in order to build a foundation for their use as quasi-topological gravitational duals. 

In this paper we obtain exact solutions for quasi-topological black holes with a U(l) electromagnetic charge in 
asymptotically anti de Sitter (AdS), de Sitter (dS), and flat spacetimes. We shall choose the parameters of quasi- 
topological gravity such that the criteria for stability and ghost-freedom are satisfied. Of interest is the effect of 
the quasi-topological parameter /i on black hole thermodynamics. After considering how the temperature-entropy 
relationship is modified, we shall examine the specific heat of the charged black holes, to see how \x allows us to 
explore a larger range of possibilities for phase-transitions than in the case of Gauss-Bonnet gravity alone. 

We find that a sufficiently negative quasi-topological parameter is able to produce uncharged spherical black hole 
solutions with two event horizons in asymptotically AdS, dS, and flat spacetimes, and that such negative corrections 
also push black hole solutions towards extremality when compared with Einstcinian solutions of the same outer horizon 
radius. We also find that the quasi-topological parameter provides greater control over the coefficients of a fifth order 
polynomial which governs phase transitions (as the denominator of the specific heat). Thus, it may be possible to 
tune the quasi-topological parameters to induce and remove phase transitions which do not occur in the Einstcinian 
case. 



II. 3 rd -ORDER QUASI-TOPOLOGICAL GRAVITY 

We consider quasi-topological additions corresponding to 3 rd -order curvature corrections to Einstein-Maxwell- 
Gauss-Bonnet gravity that maintain second-order field equations with respect to the metric under conditions of 
spherical symmetry. We use the action 



I = J <t°z^i (-2A + £, + {D _^ D _ 4) C + ^ 3 - ii^j (1) 

where D is the number of dimensions, F^ v = d^A^ , /i and A are the correction terms' coefficients, L\ = R is the Ricci 
scalar, £ 2 = R^v-ysR^ 15 ~ 4Rp V R liV + R 2 is the Gauss-Bonnet Lagrangian, and £ 3 is the 3 rd order quasi-topological 
gravity term. This term has the form 
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and is only effective in dimensions greater than D > 4, becoming trivial in 6 dimensions [15 . 
The metric ansatz we employ is 

ds 2 = -^f(r)dt* + ^f + ^dnl (3) 
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where d£l\ is the metric of a constant curvature hypersurface 
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The parameter fc = — 1, 0, 1, corresponding to hyperbolic, flat, and spherical geometries, respectively. We parametrize 
the Maxwell gauge field with a function h(r) 



A = qjh{r) 



and additionally introduce the notation j(r) = an ^ > for convenience 



_ dh(r) 
dr 

The field equations that follow from the action are 
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3 5), where k = (^g — ^prkj . These equations simplify considerably if we substitute f(r) = N 2 (r)g(r), yielding 
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III. BLACK HOLE SOLUTIONS 

We consider here the exact solutions to the field equations in five dimensions with the addition of the quasi- 
topological term. It is clear that N = 1 solves one of the field equations. We first begin with the h = case. 
Setting h = 0, it is clear that we have one equation to solve: 
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which is cubic in k. There are three branches to the solution, only one of which reduces to the Reissner-Nordstrom 
solution when X = fi = 0. The principal solution (which does not necessarily reduce to the Reissner-Nordstrom 
solution) is 
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and M is a constant of integration. This solution matches the form recently obtained for the 7-dimensional 3 rd order 
Lovelock case |16) , since the field equations are the same apart from powers of r that arise due to this difference in 
dimensionality. 



The charged solutions are obtained by solving ( 11 ) 



A = Q 



dt 
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in turn yielding 
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where Q and M are constants of integration and A is the cosmological constant. This is the 3 rd -order quasi-topological 
generalization of the Reissner-Nordstrom solution. 



The principal branch of the solution is the same as that of (13), but with 
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To ensure that the principal solution of the cubic is real, we find that the inequality 
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must be satisfied. 



Turning to look at all branches, we rewrite eq. (12 1 for g(r) as 
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where e = — AL 2 /6 is respectively equal to 1, —1,0 for asymptotically AdS, dS, and fiat solutions. The Kretschmann 
scalar diverges if the condition 



d r n — — W[k] 
dr 



dW[n] 
8k 



is satisfied. The denominator corresponds to extremal values of W[k] with respect to n. In the charged case, the 
numerator can also equal zero for 



f = 3^ S M. 



where r ns is the radius corresponding to the solution of the extremal value of W[n\. 

For simplicity we do not consider the case when the numerator reaches zero (this could result in a 0/0 solution), 
and broadly disregard solutions where W[k] reaches extremal values. Taking the derivative of equation (19), we find 
that extremal values occur when 
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To ensure that this never happens, we can force the discriminant to be negative (given that g(r) is always real as per 
restrictions above) yielding 
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We therefore see that by satisfying equation ( 18 ) we obtain meaningful roots. The value of A is all that remains to 
be chosen. For simplicity, we choose between three values of A; one for AdS, one for dS, and one for flat-space black 
hole solutions. 

In what follows, the values of A and /i shall be chosen to lie within the region where our theory has stable AdS vacua 
without ghost gravitons, and where nonsingular planar black hole solutions exist |15j . Unless otherwise stated, we 
shall use A = 0.04 and /i = —0.001. Note that for A = 0.04 we remain within this region provided /i ^ —0.000408. The 
solution ( 13|17 l is the only real solution within this region of A and fi. We note that as A, \i — > 0, this solution does 
not approach the Einsteinian black hole. Taking // — > 0, as discussed above, will lead to imaginary components in the 
function g(r) and the solution is deemed unphysical. Indeed, examining the roots of the equation ( 16 ), one can observe 
that the real root travels to ±oo as the parameters A, (i — > (depending on the direction from which is approached). 
Nonetheless, though the branch we select does not ultimately converge to the Einsteinian branch in this limit, we shall 
see that the thermodynamic behaviour does become very similar for \x = —0.001 and a small A, and so we can make 
reasonable conclusions about the effects of a quasi-topological parameter on black hole thermodynamics. In addition, 
the other two branches of the field equations only produce nonsingular solutions in regions where the AdS vacuum may 
be ghosty |T5]. Though nakedly singular solutions have been shown to have some observable characteristics [T7], we 
will choose to study black hole thermodynamics in this paper because they are more well-understood as conjectured 
duals in the gauge/gravity correspondence. 



Asymptotically anti de Sitter solutions 



We constrain A for these solutions such that they admit the limit lim r _ i . 00 g(r) — 1, using the field equation (16): 
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This restriction is introduced simply because it is convenient to compare metric functions that converge to the same 
value; the effects of the higher curvature terms can be more easily seen. When we subsequently discuss thermody- 
namics, we shall see that this definition of the cosmological constant will artificially affect the temperature/entropy 
relationship, and so it will become more convenient to compare solutions with the same cosmological constant. 

We are then able to obtain Reissner-Nordstrom solutions for a given black hole radius. These solutions are asymp- 
totically anti de Sitter provided the cosmological constant stays negative, that is, fi + A < 1. In order to meaningfully 
compare results for various cases, we consider black holes whose event horizons have the same radius, which in turn 
implies that they have different mass. 

Figure [I] shows three solutions with an event horizon of r + = 1.0 for L = 1 (L = 1 is assumed throughout the 
paper): Q = (solid), Q = 12 (dashed), and Q = 6 (dotted). It turns out that Q = 6 is near-extremal (the inset 
of Figure [I] indicates 2 horizons) , and so Reissner-Nordstrom solutions with significantly greater charge than this do 
not exist for the given event horizon radius. For Q = 12 we plot an extremal solution since we cannot obtain a 
r+ = 1.0 black hole for this case. The respective masses for Q = 0,6, 12 are M = 2.883000024, M = 7.383000024 and 
M = 18.483; just beyond the near-extremal Q = 6 case graphical methods yield an approximate extremal mass of 
M 7.3315. 



For given values of Q and one of the roots of g(r), the value of M can be obtained from eq ( 16 1. To ensure that the 
root of g(r) is the outer horizon r+, we used an iterative method to determine the given M parameter values which 
corresponded to particular roots at r = r + . To ensure accuracy in locating the roots, we required that the iteration 
algorithm not cease until convergence to within 10 -10 of our desired r + was attained. For black holes without a 
quasi-topological parameter (/i — 0), we used analytic solutions that correspond to the value of M for any given 
rightmost root r + , as previously done in the chargeless case [TS]. These solutions are listed in Appendix |A} 

Figures [2]^4] respectively illustrate solutions for k = — 1, 0, 1 of the metric functions of black holes with outer radius 
r + = 2.0, for Q values of 0, 6, and 12. We find that the functions f(r) = g(r) possess two linear zeroes (and so the 
metric has two event horizons) in the charged cases, and that they properly asymptote to 1 (though this is difficult 
to see due to the scale of the graph). 

The inlay plots show that the quasi-topological solution (with A = 0.04 and fj, = —0.001) modifies the behaviour in 
a counterintuitive manner - it does not become "more extremal" because the inner horizon does not converge towards 
the event horizon; instead for k = 1 the inner horizon actually moves to a smaller radial value, while for k = and 
k = — 1 the inner horizon does not appear to change significantly from the Einsteinian case. The higher curvature 
corrections also act to decrease the magnitude of the negative spike in g(r). 

Interestingly, for k = 1 the quasi-topological solution with Q = induces two event horizons (which can just be 
seen in Figure [4| . In this sense quasi-topological gravity induces charge-like effects into the space-time structure, a 
phenomenon that we shall see is also present in asymptotically dS and flat spacetimes. 
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FIG. 1: Comparison of Reissner-Nordstrom and uncharged Black Hole Solutions' metric functions versus radial coordinate 
for k = 0. The solutions are Q = (solid), Q = 6 (dotted), and Q = 12 (dashed). The inset is zoomed in to enhance the 
near-horizon behaviour of q(r). 
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FIG. 2: Reissner-Nordstrom Black Hole Solutions for k — —1, r+ = 2.0, where Q = 0,6,12 are in red, green, and blue, 
respectively. The dotted solutions are quasi-topological, and the solid solutions are Einsteinian. 



2. Asymptotically de Sitter solutions 



Requiring the metric functions g{r) to tend to —1 as r — ¥ oo yields from (16) 

6(l-/x + A) 



A = 
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and we shall adopt this value of A for this section. 

We can see that for A — /i > —1, the spacetime is de Sitter. The function g(r) is plotted in Figure [5] for a mass 
parameter of 0.7. In this case, the differently charged black holes do have different event horizon radii, though they 



7 




10 



FIG. 3: Reissner-Nordstrom Black Hole Solutions for k = 0, r+ = 2.0, where Q = 0, 6, 12 are in red, green, and blue, respectively. 
The dotted solutions are quasi-topological, and the solid solutions are Einsteinian. 
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FIG. 4: Reissner-Nordstrom Black Hole Solutions for k = 1, r+ = 2.0, where Q = 0, 6, 12 are in red, green, and blue, respectively. 
The dotted solutions are quasi-topological, and the solid solutions are Einsteinian. 



do not differ greatly; the event horizons are all between r+ = 0.8 and r + = 0.9. We can see that increased charge 
results in the appearance of naked singularities, and zero charge causes a negative asymptote (i.e. the solution then 
has only one horizon). 

Note that with our requirement that g(r) tends to — 1, we also need the presence of a cosmological horizon at 
some r — rc > r+, which is present in the solution of equation (13) only if k = 1; otherwise this horizon will not 
be present. In this case, the negative quasi-topological parameter modifies the solutions in the same direction as the 
positive Gauss-Bonnet parameter. These parameters also act oppositely to the AdS case - they now do induce "more 
extremal" behaviour in solutions. Sub-extremal solutions in Einstein gravity move beyond extremality for even the 
small values of A and fi we consider; specifically, Q — 0.6 is not extremal in the Einsteinian case, but is seen from fig. 
[5] to be nakedly singular for both Gauss-Bonnet and quasi-topological solutions, where only one horizon exists. 



Recall that for a value of A = 0.04, the condition on [i producing a ghost-free AdS vacuum is [i < —0.000408. 
We find in this asymptotically dS context that for values of (i more negative than this we obtain solutions with 
three horizons even for Q = 0. In this manner the gravitational effects of charge thus are induced by the higher 
curvature corrections in quasi-topological gravity. However a sufficiently negative value of fj, yields solutions with a 
naked singularity and one cosmological horizon. For example a value of /j, m —0.011 causes both the Q = 0.2 and 
Q = solutions to move past extremality, yielding singularities surrounded by a cosmological horizon only. We note 
that in this example, /i < —0.001 with A = 0.04 does not in the asymptotically AdS context guarantee positivity of 
energy flux in the dual CFT [7J. In the asymptotically dS context we are considering in this section we find that 
violating this bound on /i can yield nakedly singular solutions. 

Attempting to produce for Q = a space-time with only 2 horizons entails increasing /i above the region where our 
solution is real. Hence the "charge" -like behaviour induced by the quasi-topological term is a result of our branch not 
reducing to the Einstcinian branch. 
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FIG. 5: Solutions in de Sitter space for k = 1, where Q = 0,0.2,0.6,0.7 are in black, green, red, and blue, respectively. Here 
A = 0.04 and [i = —0.001. Quasi-topological gravity is dotted, Gauss-Bonnet is dashed, and Einsteinian is solid. 



3. Asymptotically flat solutions 

Removing the cosmological constant (setting A = 0) will result in g(r) being asymptotically cx L 2 /r 2 , rendering the 
spacetime asymptotically flat. Again, we also require that k = 1 in order to obtain asymptotically Minkowski space 
with our solution |l3]). We have plotted our solution in Figure [6 as a class of solutions of a mass of M — 0.8. 

Similar to the dS situation, we notice the return to a two-horizon solution when a negative quasi-topological 
parameter is introduced, without charge. We also notice that for Q = 0.6, the solution is a naked singularity for 
Einsteinian, Gauss-Bonnet, and quasi-topological solutions, and the negative quasi-topological term acts in the same 
manner as a positive Gauss-Bonnet term. 



IV. RN THERMODYNAMICS 



We compute the entropy and temperature of the Reissner-Nordstrom solutions in order to determine their stability. 
In addition to examining AdS solutions (where A = — 6/L 2 ), we have also performed the analysis for asymptotically 
de Sitter and asymptotically flat solutions, with A = 6/L 2 and A = 0, respectively. The reason we now change our 
convention and examine a cosmological constant such that the metric functions g(r) do not converge to the same 
value in quasi-topological, Gauss-Bonnet, and Einsteinian gravity is because the thermodynamics are sensitive to the 
value of the cosmological constant that the metric sees, and we want to study the thermodynamical behaviour of the 
black hole solutions under the effect of a single cosmological constant. 
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FIG. 6: Solutions in asymptotically flat space for k = 1, where Q = 0, 0.2, 0.6, 0.7 are in black, green, red, and blue, respectively. 
Here A = 0.04 and fj, — —0.001. Quasi-topological gravity is dotted, Gauss-Bonnet is dashed, and Einsteinian is solid. 



We compute the entropy from the Iyer/Wald prescription |19j where we assume the ansatz 



A A . L 2 . , 2 L 4 



Sk = Jr U + 6A/C— +9 M fc l —\ (21) 

where D is the number of dimensions and A is the surface area of the black hole. We use this entropy to check the 
second law of thermodynamics. The temperature is determined by enforcing regularity of the Euclidean section of 
the spacetimc; 
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r =UH„, (22) 

where is the radial coordinate of the outermost horizon of the black hole. 

Turning next to the first law of thermodynamics for charged quasi-topological black holes, recall that the first law 
supplies the equations 

w = T « 23 » 

<J>Jk Q=const 
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where M signifies the mass parameter evaluated at horizon (all quantities are at horizon here). We have verified that 
this law holds for all of our solutions. We note that <I> is the total electric flux over the horizon, evaluated by 

s 



$ = / (TxA r (25) 



Since curvature doesn't affect the global electrostatic potential (as it's measured at infinity with respect to the horizon), 
we are able to use the same result as [20]. That is, 

* = 16 ^ (26) 



To compute these values we can differentiate the entire expression (16) and solve for g'(r). Substituting values for 
the mass parameter and charge parameter at the horizon, we are able to absorb leading coefficients into the values 
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for mass and charge, giving us "true mass" and "true charge" m, q. We find 

M 
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A negative slope in a temperature-entropy plot produces a situation in which the black hole will be thermody- 
namically unstable as its temperature increases with decreasing radius (decreasing entropy), as seen previously for 
Einsteinian AdS [21]. In the AdS case this causes small black holes to either decay via Hawking radiation or else 
tunnel into larger black holes until thermodynamic equilibrium with the background is reached. However for AdS 
black holes of sufficiently large charge, we see from Figure[7]that no phase transition takes place as the slope is positive 
in all cases considered. Similar to the effect of the quasi-topological term in numerical black hole solutions [T2], the 
thermodynamic behaviour of the quasi-topological term closely follows Gauss-Bonnet thermodynamics for black holes 
beyond a certain radius (on the order of ~ 1). 
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FIG. 7: log (Temperature) versus log(Entropy) for Quasi-topological, Gauss-Bonnet, and Einsteinian (dotted, dashed, solid) 
asymptotically AdS RN Black Holes of k — — 1, 0, 1 (red, blue, green) and Q = 6, A = 0.04, fi — —0.001; note that the dotted 
and dashed lines yield almost identical curves for k = — 1 and that all three curves are effectively indistinguishable for k = 0. 



In the limit of small black holes, upon closer inspection we see a deviation between the quasi-topological solution 
and the Gauss-Bonnet solution. We terminate our plots in the limit of extremality for the black holes. For [i — —0.001, 
it is apparent that the quasi-topological term does not induce much change over the Gauss-Bonnet case. We find that 
a larger negative value of \i induces more dramatic behaviour, where for k = and k = 1 the slope of the curve in 
a logT vs. logs' plot asymptotes to horizontal in the limit of small black holes. While this behaviour is illustrative 
of the effect of [i, it must be taken cautiously because large negative fj, no longer guarantees positivity of the energy 
flux in the dual CFT [7]. 

We have also plotted solutions for the asymptotically fiat and asymptotically de Sitter (Figure [8]) cases. As 
previously mentioned, only the k = 1 cases are physical. Because the de Sitter black holes exist dynamically with 
a cosmological horizon, an outer horizon, and an inner horizon, we choose to measure black hole temperature on 
the outer horizon. In addition, as the black hole radius of the de Sitter solutions becomes large, the outer horizon 
approaches the cosmological horizon (forming a Nariai solution [52]), beyond which the space-time has a singularity 
with only a cosmological horizon. The upper bound for the dS case in Figure [H] is the Nariai limit. 

In the case of asymptotically flat black holes, we see that our results exhibit close agreement with ref. [23] . 
There, the same behaviour regarding uncharged higher-curvature solutions was found - not unexpectedly, as the field 
equations are similar. In addition, their thermodynamic analysis showed that the black holes have a region of stability 
for r extremal < r + < r unsta bie, commensurate with our analysis. 

We see in Figure [8] that the dS black holes undergo a single phase transition from stable to unstable as their horizon 
radius grows. At this transition point (where the curve in Figure ^ reaches a maximum), the black hole outer radius 
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FIG. 8: log(Temperature) versus log(Entropy) for Quasi-topological, Gauss-Bonnet, and Einsteinian (dotted, dashed, solid) 
asymptotically de Sitter (blue) and asymptotically flat (red) RN Black Holes of k = 1; Q = 0.6, A = 0.04, /i = —0.001. 
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FIG. 9: log(Temperature) versus log(Entropy) for asymptotically anti de Sitter Reissner-Nordstrom Black Holes of Q — 6 
(dotted), Q = 3 (dashed), Q = 0.5 (dash-dot), Q — 0.1 (long-dash), and Q = (solid) for k = —1,0,1 (red, blue, green); 
A = 0.04, fj, = —0.001. 

is r + sa 0.485 for the Gauss-Bonnet case. For the flat case, another phase transition is seen in Figure [8] for k = 1, 
occurring at the black hole outer radius r + s» 0.580 for the Gauss-Bonnet case. In the quasi-topological case, these 
numbers change to r + m 0.496 and r + w 0.600, respectively. 

In Figure |9j we illustrate the effect of the charge parameter Q on the quasi-topological gravity solutions. The solid 
lines are uncharged black holes (and exhibit thermodynamic behaviour described in |12j). We see that increasing the 
charge produces much more stable black holes with larger radii, where a very large increase in temperature is required 
to deliver a small increase in entropy. This behaviour is precisely what is expected when increasing the charge - 
extremal black holes will have larger mass. Furthermore, as Q increases the distinction between the different values 
of k becomes less pronounced, being nearly indistinguishable for Q = 6. 
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4- Specific Heat 

The specific heat is given by 



Computing this quantity 



ft = f (*>) 



ft - < 30 > 

or + or + 



we find that roots exist, and therefore there are unstable regions. Recall that a negative specific heat will result in a 
thermodynamic instability with the vacuum |24| . For k = 0, the only real roots of the specific heat are 



r h =±l/2 V V ^ (31) 

Therefore there are no unstable regions unless A is negative, as is required in order to obtain valid black hole solutions 

for k — 0; when A = —6 (AdS), we find that 77, = ® &1/B is a transition. For Q — 6, this gives us a phase transition at 
0.953. 

For k = 1, the equation for specific heat becomes significantly more complex, and characterizing solutions is much 
more difficult. We begin by looking at the case for A = —6, k = 1. Direct substitution leads to the expression 

12tt (-48 r 6 - 24 r 4 - 24 fi + Q 2 ) (-r 4 - 2 A r 2 + 3 fj) 2 

v = r (-288 r 8 A + 720 r 6 fi + 144 Ji X r 2 - 6 Q 2 X r 2 - 48 r 6 A + 336 r 4 /i + 3 Q 2 fi - 48 r 1Q + 24 r 8 - 72 /j 2 - 5 Q 2 r 4 ) 

(••52) 

From the numerator, the only non-saddle zeroes will result from the vanishing of the term (Q 2 — 48r^ — 24/i — 24r 4 J . 
We see that a transition therefore occurs under the condition 



Q = ±2^/127-6 + 6r 4 +6/j (33) 

where r/j is the radius of the transition, denoting r + by when it corresponds a phase transition radius. 

The denominator produces a tenth-order polynomial whose roots correspond to divergences of the specific heat, 
resulting in complicated phase transition behaviour. This type of divergence also occurs in, for example, studies of 
quantum gravity |25j . While explicit (numerical) computation of these roots is a challenging task, we can see that 
roots exist for the value of Q such that 



Q = ±2 ^V (SM-GA^-Sr, 4 ) ( 34 ) 

We need to make sure that this corresponds to a realistic (i.e. real) value of Q. For the values of [i < and A > 
in the stable ghost-free region that we consider, the denominator under the square root is negative and so we must 
find values for rt such that the numerator under the square root is also negative. It is clear that for sufficiently small 
|A| and that this is always possible - for example A = 0.015 and /j, = —0.0005 there exist real solutions for Q. 
However for A = 0.04 and fi = —0.001 the numerator never becomes positive, and therefore, there are no divergences 
in the specific heat. 

As another example, consider a cosmological constant (A = —1) for which r/j = 0.8031577349. This results in 
Q = 1, a valid solution. The specific heat is positive and increasing at r + < r^, and the black hole undergoes a 
phase transition where the specific heat becomes singular at r + — r^. As r + > r^, the specific heat increases (though 
it stays negative), until r + w 0.92 whereupon the specific heat begins to decrease again. Another root can also be 
found: r/j = 1.59885311 also results in Q = 1, and this is identified with another transition, where the negative specific 
heat diverges to —00 and reappears at r + > 1.59885311 with positive value. As Q increases beyond a critical point 
(somewhere between Q = 2 and Q = 3), the discontinuities in the specific heat disappear (there are no longer any 
real positive roots in the denominator), and we have solutions similar to the A = —6 case. 

As a verification that our thermodynamics plots are correct, we have checked that the maxima of our curves 
correspond to divergences in the specific heat. For example, for Q = 6 there is a transition from negative specific heat 
to positive specific heat at r + = 0.877 in the AdS quasi-topological case (k = 1). Above that radius, the specific heat 
stays positive, so our specific heat results agree with the thermodynamics in Figure [7] This transition is unphysical - 
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it happens to be the extremal radius for the black hole and, not coincident ally, it is also the radius at which we stop 
our temperature-entropy plot of Figure [7J 

Switching to A = 6 for de Sitter solutions, the function for specific heat is similar, and we solve for the existence 
of roots for fi = —0.001 and A = 0.04 in exactly the same manner as for AdS solutions. At Q — 0.6, we find that 
the specific heat crosses from negative to positive for these values of (/x, A) at r + fx 0.39, followed by a divergence 
to oo at r + ss 0.4945. This implies three phase transitions. The first and third of these 'transitions' are where 
the outer horizon respectively merges with the inner and cosmological horizons; both are left off of the plot. The 
meaningful phase transition is the one for which the black hole transits from instability (where the temperature (and 
size) decreases with an increase in entropy) and a stable phase, where the black hole's temperature increases with an 
increase in entropy. The single phase transition can be seen in Figure [S] where the curve reaches a maximum. Similar 
behaviour occurs in the Gauss-Bonnet and Einstein cases for correspondingly smaller values of r + , as is clear from 
Figure [8] This type of discontinuity in the specific heat has also been studied for Born-Infeld electrodynamics [35] . 

Finally, flat space (A = 0) solutions also lead to similar occurrences. For Q = 0.6, /i = —0.001, and A = 0.04, the 
specific heat will cross from negative to positive around r s» 0.35, followed by a divergence to oo at around r s» 0.59. 
The specific heat will then increase from — oo, eventually reaching a turning point and decreasing (never crossing 
zero) . This implies two phase transitions. Again, the first transition is where the black hole becomes extremal and is 
left off of the plot. The remaining phase transition is again the black hole's transition from an unstable phase (where 
they possess large horizon radii) to a stable phase (for black holes with small horizon radii) . This is again seen in our 
temperature/entropy plot of Figure [7j 

Summarizing, the intriguing result is that we effectively have control of a 5 th order polynomial with coefficients 
that depend on A, /i, A, and Q, allowing us to theoretically (without constraints on /x or A) produce up to five 
discontinuities in the specific heat. In general, this polynomial will contain five important coefficients and one scale 
factor (which is unimportant when computing its zeroes). In the Gauss-Bonnet scenario we have three parameters 
upon which four coefficients in the polynomial depend, before A is fixed, as the constant term disappears (group the 
denominator of the specific heat into terms of r 2 , r 4 , r 6 , etc.). For the quasi-topological case we have four parameters 
to fix five coefficients. An additional quasi-topological term [B] might allow us to obtain the final parameter, and 
perhaps it will form an independent basis of these parameters on each term, giving us access to a maximum of five 
phase transitions (that is, allowing us to redefine four unique parameters eti, aa, «3, a-4, for which the denominator 
of the specific heat appears as r 10 + air 8 + a^r 6 + a^r 4 + a^r 2 + as). 

While adding a fourth-order quasi-topological term appears to be a useful endeavour, we leave its effect open to 
question since it is not immediately clear, if successful, what utility a complete basis for the roots would be, given that 
a ) the solution to the quintic is not in general expressible in radicals and b ) the parameters fi, A, Q all have various 
restrictions and so they do not extend over the reals. In addition, the possibility exists that a quartic quasi-topological 
term will increase the order of the polynomial (though we do not expect this to be the case, as it did not increase 
from Gauss-Bonnet to third order quasi-topological gravity). 



V. CHARGED BLACK HOLES IN iY-TH ORDER QUASI-TOPOLOGICAL GRAVITY 

In this section we consider K-th order quasi-topological gravity in (AT + 1) dimensions. The specific form of the 
Lagrangian for K-th order quasi-topological gravity has only been obtained for K = 3 [15] and K = 4 [27] , However 
we can make some remarks for general K based on this work, even though the specific Lagrangian has yet not been 
found. 

Consider the spherically symmetric ansatz 

,y>2 7~2 J™2 „2 

ds 2 = - T2 N\r)g{r)dt 2 + + -^l (35) 

A = h(r)dt (36) 

where dT, 2 , is the metric of a (N — l)-dimensional constant curvature hypersurfacc. 
Maxwell's equations yield 



r N-l 

N 



W') =0 (37) 



implying 



A = (N-2)r*- 2 dt (38) 
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as the solution for the gauge field. Since the field equations can only be 2 nd order in the metric functions in the 
spherically symmetric case, it is reasonable to conjecture that the field equation for the metric function g(r) is 



A' 



,N 



„N-1 



fe=0 



Q 



,N-1 



(39) 



for K-th order quasi-topological gravity in (N + 1) dimensions, where K = (kL 2 r 2 — g). We easily obtain 

A 



T 1.9] = 



M 



i=0 



r N (N - 2)r 2N - 2 

and the metric function for Einsteinian gravity (setting fii = 1) is 

M 



9 2 



■k- 



Q 



L 2 r N-2 L 2 (N-2)r 2N ~ i 



(40) 



(41) 



This equation has the same form as the corresponding Q = situation in Lovelock gravity [28l [29], where K < [^] . 
In quasi-topological gravity K is not restricted (though there are certain values of N for which the construction does 
not make sense). 

We start by considering the case where k = 0. We find that the planar solution of equation (41 ) admits a vanishing 
g(r) only when 



Mo 



M 

N 



Q 2 

(N - 2) r h 2N ' 2 







or corresponding to roots of 



Mr h N ~ 2 



2N-2 



Q 2 



N -2 



(42) 



(43) 



These roots can be numerically obtained. We can see that in higher dimensions, the high polynomial degree results 
in a rich dependency of horizon solutions on M and Q. 
Now, more generally (for any k), at the horizon we have 



kl7 

Th 2 



so that 



Mr h N - 2 



fJ-orh 



2N-2 



Q 2 



N-2 



(44) 



(45) 



This is not analytically solvable but the branch we are interested in should be a perturbation of the branch for the 
uncharged case |28j . Our solution will be a continuous transform of this case and so we also relate our analysis to the 
Einstein-Hilbert (EH) branch. 

For a sufficiently small value of Q, our Gauss-Bonnet EH- branch is 



1 

2A 



2XkL 2 



-W1-4A 1- 



M 



3r 4 I" 24 r 6 



(46) 



We can see that this is equivalent to our previously obtained solution for g(r) as given in Appendix [A] This also 
happens to be the branch of the solution connected to the Einsteinian solution; as A — » 0, we recover the Einsteinian 
solution (A3). 



VI. DISCUSSION 



In this paper we have examined Reissner-Nordstrom-type black holes in quasi-topological gravity. We have intro- 
duced third-order quasi-topological curvature terms to the Lagrangian, and we find that the introduction of a Maxwell 
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field leads to AdS solutions which have the double-horizon characteristic of the Reissner-Nordstrom black hole. As 
expected, the small quasi-topological terms do not drastically alter the solutions. 

We have seen that not only can the Maxwell charge parameter induce two horizons in flat and de Sitter spacetimes 
(as expected), but a quasi-topological parameter - within the stable ghost-free region of parameter space - is able to 
create two horizons without any charge parameter present. In addition, along with the Gauss-Bonnet parameter, the 
higher curvature correction terms are able to take a near-extremal black hole in Einstein gravity into a singularity 
with no horizons, in flat space. 

We have also provided further verification that a small, negative quasi-topological term will act similarly to the ef- 
fect of a positive Gauss-Bonnet term with respect to stability. We have also seen how the quasi-topological parameter 
can affect phase transitions. Because the quasi-topological parameter is restricted to small values, we are not able 
to add or remove phase transitions at will, but we do see that for specific cases, fj, is able to broaden the region of 
thcrmodynamical parameter space that we are able to probe. Specifically, the specific heat divergences are governed 
by a quintic polynomial in black hole radius. The Gauss-Bonnet solution allows us to control two coefficients inde- 
pendently, while the quasi-topological term gives us control over three of the five coefficients (once the cosmological 
constant had been determined). 

Regarding future work, Tarrio and Vandoren |30| have shown parallels to Reissner-Nordstrom black holes in gravity 
with Lifshitz scaling symmetries, with unexpected behaviour. It would be interesting to see the modifications to these 
behaviours in a Lifshitz symmetric gravitational theory |31[ I32j with higher order quasi-topological terms. Another 
interesting phenomenon is the surprising result of instability of RN black holes in dimensions greater than 6 |33] . Will 
this result hold in the quasi-topological regime? 

Appendix A: Analytic Asymptotically Anti de Sitter Solutions for M given r+ 

For Gauss-Bonnet gravity, we obtain a simplified solution for f(r) — g(r) of 



, , 6r 3 + !2XkL 2 r - J36r 6 - 144Ar 6 - 6XL 8 Q 2 + 48XMr 2 , , , 

9( r ) = ( A1 ) 



12Ar 3 

and, setting L = 1 and solving roots for M, we see that 



-24r+ 6 - 24Ar+ 6 + Q 2 + 24r+ 4 fc + 24Xk 2 r + 2 
M = — ~ (A2) 



Similarly in Einsteinian gravity, 



with roots 



, , L 8 Q 2 - 8Mr 2 + 24r 6 + 24kL 2 r 4 , s 

9(r) = ^6 ( A3 ) 



M = Q2 + 24 -+ 6 + 24 ^ (A4) 
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